Uncertainty Principle Enhanced Pairing Correlations in 
Projected Fermi Systems Near Half Filling 

B. Sriram Shastry 
Indian Institute of Science, Bangalore 560012, INDIA 

Abstract 

We point out the curious phenomenon of order by projection in a class of lat- 
tice Fermi systems near half filling. Enhanced pairing correlations of extended 
s-wave Cooper pairs result from the process of projecting out s-wave Cooper 
pairs, with negligible effect on the ground state energy. The Hubbard model is 
a particularly nice example of the above phenomenon, which is revealed with 
the use of rigorous inequalities including the Uncertainty Principle Inequality 
In addition, we present numerical evidence that at half filling, a related but 
simplified model shows ODLRO of extended s-wave Cooper pairs. 
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There is considerable current interest in the possibility of purely electronic interactions 
driven superconductivity as a mechanism to explain the High T c superconductors. While 
it is well known for the uniform electron gas that purely Coulomb repulsion terms lead to 
superconductivity in higher angular momentum channels JIJ, albeit with very low transition 
temperatures, here the search, guided by experiments, is predominantly for single band mod- 
els that display such behavior in the proximity of half filling on a lattice. The prototypical 
example is that of the Hubbard model [|J |§, although its tendency (or otherwise) towards 
superconductivity remains an unsettled issue. 

In this work, we study a class of many body Fermi systems on a lattice, under the 
influence of a projection of s-wave Cooper pairs. Recall that one has an inhibition of s-wave 
ordering within weak coupling BCS theory for models for on-site repulsion in addition to 
the usual phononic coupling. In contrast, we project out s-wave Cooper pairs in the present 
work. The study of most projected models, generally justified by their status as "fixed 
point" Hamiltonians in some underlying scaling theory, has been a rich source of new and 
interesting models in the field of correlated fermi systems. A prototype is the Gutzwiller 
wavefunction, wherein upon removing double occupancy, effects such as enchanced effective 
masses follow near half filling, and these are crucial in our understanding of almost localized 
fermi liquids Q. At half filling we find insulating wave functions with enhanced spin-spin 
correlations ||, that are regarded as typical of Quantum Spin Systems in low dimensions 
with S — 1/2. At the level of the Hamiltonian, projection leads to interesting new models, 
such as the various limits of the Hubbard model, e.g. large U giving the t — J model, U = oo 
giving the Nagaoka limit, and several examples in single impurity models. It seems worth 
remarking that projection is a theoretical device that is genuinely strong coupling and non 
perturbative, making it difficult to treat with conventional methods. In the present work 
the consequences of s-wave projection are found by combining a set of known inequalities in 
a novel fashion, and lead to surprising insights detailed below. 

We consider the model defined on a d-dimensional hyper cubic lattice with a Hamiltonian 
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H = T + UY,n i ^n li + U s B^B (1) 

i 

where T is the kinetic energy term Ylka e (k)c\(k)c a (k) . The second term is the Hubbard 
repulsion term (i.e. U > 0). We will consider other forms of interaction below, but the 
argument is simplest for the Hubbard interaction written above. We will take the number of 
sites as £ and denote the density of particles by p = N/C We will also denote H = H — fiN, 
where p is the chemical potential and N the number operator. The third term is new, with 
the operator B = J2i exp{i0j} bj and bj = Cj^Cj^. If we take U s to be 0(1/ £), and negative 
then this term would, in a weak coupling BCS like theory, inhibit the formation of s-wave 
Cooper pairs. The coupling constant U s is taken of 0(1) and positive in the present work 
and corresponds to projecting out the appropriate Cooper pairs, at general fillings. Precisely 
at half filling, the influence of the new term is more subtle as noted later in the paper. 
Although various choices of the phase angle (pi generate different examples, two of the 
interesting ones are (a) (pi = leads to a suppression of pure s-wave superconductivity, 
and (P) (pi = r*j ■ {n,ir, ..} suppresses the so called eta pairing 0. We will also consider 
a third possibility (7) obtained by setting B = Sfc ex P{^0(^)} b(k) where b(k) = c^^c^ 
and with an arbitrary function (p(k) which can be used to vary the relative phases between 
different momenta. This last class of operators, however forces us to the case of U = 0, in 
order to obtain any results. The case (a) appears to be the most interesting physically, 
the others are included for completeness. 

We first note that for lattices that are bipartite, and where the electronic hopping only 
connects unlike sublattices, we can make a particle hole transformation Cjj — > ( — l) 6i cL, with 
9i = 0, 1 for the two sublattices, whereby the energy satisfies E[U, U s , p] = E[U,U S ,1 — p] — 
C(l — p)(U + U s ). At half filling (p = 1) the chemical potentials for adding and subtracting 
a particle add up as: p + + \i_ = U + U s . At this filling, the new term U s plays a crucial 
role in allowing doubly occupied sites and holes to wander away from each other, and infact 
encourages charge fluctuations, whereby the usual Mott insulating state of the Hubbard 
model at half filling is heavily discouraged. 
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We now use a simple but useful inequality |7j 

<^fj \M^ [H,M]\^ >> (2) 

where > is the ground state of H, and M is an arbitrary operator. Using M = B we 
find on using the important commutator [B,B^\ = C — N, valid in all cases (a), ((3) and 
(7), that 

< B ] A >> {U s (C-N + 2) -2(i + U} < B ] B > . (3) 

In the case of (7), the above Inequality holds only with U — 0. Note that the LHS of above 
is forced to be real and to be positive from the Inequality We denote ground state averages 
by angular brackets as above, and the operator A is given by A = [T, B] . For the two cases 
of the phase fa in Eq© , (a) A = -2 £ fc e(k)b(k), and (/?) A = - £{e(£) + e(k + II)} b(k) 
with II = {ir, 7T, ..}. In the popular case of nearest neighbour hopping on the hypercubic 
lattice, (a) corresponds to the extended s-wave pairing operator, whereas gives zero. 
A non zero result is obtained in the latter case only when the hopping connects sites on 
the same sublattice. In two dimensions, for example, with e(k) = —2t(cos(k x ) + cos(k y )) — 
2t' cos(k x ) cos(ky), we find A = 4t' J2( cos (k x ) cos(k y )) b(k). 

We now use the Cauchy-Schwartz inequality to bound the LHS of Inequality(||) as 

< A^A >< B ] B >>< B^A > 2 . (4) 

Combining Ineqs(|| ,f| ) we find 

< A* A > > {U S (C - iV + 2) - 2/1 + U} 2 < B^B > . (5) 

Again note that in case of (7), Ineq(|5]) is valid only with U — 0. We note that in the RHS 
of Ineq(||) the prefactor is of the 0(C 2 ) provided we are at a thermodynamic filling p < 1. 
Exactly at half filling the inequality is less useful. At any filling p < 1, we can deduce that 
< B^B > is very small. In fact we will show that it is o(C) rather than 0(C). If it were of 
the 0(C), ( as indeed it is in the ground state of the free fermi gas), then < A* A > has to 
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exceed a trivial upper bound of the 0(£ 2 ) ||. If < B^B > is of the 0(1) then we have two 
consequences that are mutually incompatible (at least when U = 0). To see this, assume 
that < B^B > is of the 0(1) and so we find from the Feynman Hellman Theorem 
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E[U, U„ p] = E[U, 0, p] + f ' dU' s < B^B >u> (6) 

J 

= E[U,0,p] + o(C). (7) 



The other consequence of Ineq(^) is that < A^A >~ 0(£ 2 ), i.e. we have Long Ranged 
Order (ODLRO) in the operator A ||. This is possible only if the energy increases by 
terms of the 0(£), at least in the case when U = as is seen from a diagonalization of a 



bilinear Hamiltonian adding the kinetic energy T and A with coefficients of the 0(1) JTU 
A consistent possibility |TT 



is 



< B ] B >= 0(l/£) and < A ] A >= 0(C), (8) 

along with Eq(|7|). One immediate consequence of this result is that the energy per site of 
the model in Eq(JI|) Lim c ^ooE[U, U s ,p}/£ is identical to that of the pure Hubbard model ( 
i.e. U s = 0) at all U or filling p ^ 1. 

Another important consequence is that the chemical potential is unchanged by U s until 
we reach half filling (/i = BE /ON) and therefore the compressibility is unchanged by U s 
(since 1/k— Np(dp:/dN)c). At precisely half filling, the chemical potential jumps and the 
compressibility vanishes. The value of /x at half filling for the case of bipartite symmetry 
was given above as (U + U s )/2. 

We see that the suppression of a correlation of the type < B^B > occurs with remarkable 
efficiency through Ineq(§). When we recall that < [5, ST] >= £(1 — p), it is seen that 
the fluctuations of B + B^ in the ground state diminish on approaching half filling, i.e. 

< (B + £> T ) 2 /£ >= 1 — p. This immediately suggests that "conjugate variables" in the sense 
of the Uncertainty Principle, should exhibit enhancements by similar factors. The following 
form of the Uncertainty Principle is most useful, for any two operators a and b ( such that 

< a 2 >= =< b 2 > and [a, b] = 0) we have |L2| 



< a)a + aa) > < Vb + bb ] >> \< [a) , b] > | 2 . (9) 

We now use this with a — > A and b —>■ B, and also the results [A, B*\ = IT and [A, A*] = 
Xa = 4E^W 2 (1 - Ea4(k)c a (k)) to find 

Both terms of the RHS of the above Inequality are of the 0(C), and the second term 
remains bounded as we approach half filling, infact vanishing for the case of a symmetric 
band around zero energy. This implies that the first term dominates and hence we conclude 
that the correlation function < A'A > grows without limit as half filling is approached. We 
should remark that any operator of the type [T, [T, ..[T, B}..] in the place of A would end up 
having similar enhancements in its correlations, since it would be a bilinear in the c* and 
have similar commutation with B 



We next consider other kinds of interactions, different from the Hubbard model [ I3| . In 
this case, we can still use lneq(0) and also Ineq(f|) to find in place of Ineq([|), 

< F ] F >> {U s (£-N + 2) -2/j} 2 <B^B>, (11) 

where F = A + C with C = -[B, V int ], so that F = [H, B}. The norm on the LHS of Ineq(|TT]) 



can be bounded by the triangle inequality as < F^F >< [V < A^A > + V< CtC >] 2 and 
hence we need, in addition to the previous estimates, one of < C^C >. This of course 
depends upon the nature of the two particle interaction, and has to be examined for each 
model separately. However, for "generic" repulsive short ranged models, it seems clear that 
this object, like < A^A > should be bounded from above by a number of 0(£ 2 ). With this 
assumption, the remaining argument goes exactly as in the case of the Hubbard model, and 
we again conclude that < B'B > is at least as small as o(C), and in fact probably 0(1/ C), 
and that the ground state energy is as in Eq(^). The uncertainty relation InegflTUP needs 
only the faircekly weak first condition < B^B >~ o(C), and hence we conclude that the 
mechanism of order by projection works for generic short ranged repulsive models near half 
filling. 
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We have thus found enhanced correlations as we approach half filling, and by continuity, 
we may expect ODLRO in the operator A. The inequalities given above do not constrain 
correlations sufficiently, and we turn to other methods. Before doing that, we introduce a 
simpler version of the models above, namely 

H s = E(e„ - + 1) + U S £ a+a~, (12) 

n=l n,m=l 

where a z etc are the usual Pauli matrices, and e n are an ascending set of energies. This 
model is intimately related to the U — version of our starting problem Eq([|), using the 
pseudo spin representation cr| + 1 = ^2 c n CT {kj) and = c^(kj)c}(—kj), in the subspace 
where both (k, f) and (—k, j) are simultaneously present or absent. The Hamiltonian Eq(|l]) 
commutes with the operator v = J2k n ^{k) n i{~k)i an d its operation is identical to that of H s 
provided we specialize to various sectors labeled by the eigenvalues of v (0 < < v > < N/2), 
and further choose appropriate degeneracies for the energies. We simplify by choosing our 
energies in H s above to be non degenerate, and pick them to be e n = {n— (L+l)/2}/(L— 1) so 
that the band is symmetric about zero and the bandwidth is unity. Each up spin corresponds 
to two ( fermi) particles of the original problem. The filling in this problem is clearly 
p = N/L with N = Sj(c z (i) + 1) • The chemical potential at half filling is U s /2 by particle 
hole symmetry. 

The model can also be viewed as a lattice of N/2 hard core particles sitting in a constant 
electric field that tries to localize them in regions of low potential, and an infinite ranged 
hopping that tries to delocalize them. The results proved for the starting Eq([l]) namely 
Ineqs(|5|,p!0|,|8|), are equally true in this one dimensional spin model, provided we identify 
B = J2j a j an d A = —2 J2j ■ Away from half filling, i.e. when af ot ^ 0, we see that even 
in the limit of large U s , there is a large number of states, infact states with S tot = L(l — p)/2 
and Sl ot = —L(l — p)/2 i.e. highest weight states of the rotation group, which have a 
null eigenvalue of the hopping term f/ s E f^m- The Zeeman energy term has non zero 
matrix elements within this manifold. In the case of half filling p = 1, the Zeeman term 
necessarily connects singlet states with triplets and hence the energy is unable to escape 
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the influence of U s . At half filling and for large U s , we can use degenerate perturbation 
theory to find an effective Hamiltonian to lowest order in 1/U S . To do this we consider 
the action of H s in Eq(|T2|) on the space of l Cl/ 2 / (L/2 + 1) singlets spanned, for example 
by the Non Crossing Rumer Diagrams [14]]. A typical non orthogonal state is given by 
ipp = [Pi, P2]- ■ ■ ■ [Pl-i, Pl]- where P is one of the permutations of the set {1, 2, ... , L} 
giving a Non Crossing Rumer Diagram, and = (ctiPj =F A a i)/V2 is a singlet (triplet) 

with s z = 0. The action of the operator Eqfll2f) can be projected into this subspace, by 
using the relation [l,2]+[3,4]+ = |(l-2ni 3 )[l, 2]_[3, 4]_ +^ 9 " inte * with % the permutation 
operator, and leads to the following Quantum Dimer problem: 

H qd ipp = — ^(e P . - e Pj+1 f i) P - 
s 3 

s j+Kk 

This model is quite non trivial to work with, but does reveal that the diagonal terms favour 
singlet bonds that connect the largest energy separations, and the mixing terms oblige us 
to take non trivial linear combinations in this space. 

We study the interesting half filled limit by studying the sector a^ otal = of Eq(|12"D 
directly. We diagonalized the problem numerically for chains of length up to 14, and studied 
the ground state energy as well as the correlation function < A^A >. It is clear that a non 
zero extrapolation of T to a number of the 0(1) would imply ODLRO in the A field. In the 
figure we plot the parameter r = j- < A^A > / < A^A > non for three values of U s (= 2, 4, 8). 
The data seems to be consistent with this hypothesis, and fits well to V = T 00 + \a\/L± \b\/L 2 
, with non zero T^. In the inset of the figure, the ground state energy per site is plotted 
for the same three values of U s against 1/L, showing that the energy does depend on the 
coupling at half filling, implying that the U s term cannot be viewed as a projection at 
this particular filling. The dependence is consistent with finite sized scaling with a form 
E/L = eoo + \a\/L - \b\/L 2 + 0(1/L 3 ). 

At half filling, the new model Eq(|IJ) is almost certainly non-insulating, and likely to 
be superconducting in a complementary pairing state. The presence of hopping terms for 
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the doubly occupied sites makes their number density nonzero, unlike in the pure Hubbard 
model. By continuity in filling, we expect the pairing correlations to be divergent for any U. 
Our numerical results for the reduced model, the spin model of Eq([T^) are consistent with 
ODLRO at half filling. It is not, however, straightforward to write down a mean field theory 
that captures the correct ordering in the model, since the Hamiltonian does not contain 
explicit terms that favour any kind of ordering, these are generated by the dynamics rather 
indirectly. 

In summary we have seen that the effect of projecting out s-wave Cooper pairs in a class 
of Fermi systems leads to surprising results. The ground state of the projected model may be 
viewed as being essentially degenerate with that of the original model and yet the extended 
s-wave pairing correlations are hugely enhanced near half filling. This effect, namely order 
by projection requires a lattice Fermi system near half filling to occur, and has no natural 
counterpart in continuum Fermi systems. In this regard, as well as in the form of the 
enhancements 1/(1 — p), it resembles the results of the almost localized fermi systems 
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Figure Captions 

• Long Ranged Order parameter r = t< A^A >/< A^A > non versus 1/L for U s = 2, 4, 8 
(bottom to top) chains of length 4, 6, .., 14 at p = 1. The inset shows the Ground State 
Energy per site for the same values of U s versus 1/L. (bottom to top). 
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